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GLOBAL ATTRACTOR FOR A GINZBURG-LANDAU TYPE 
MODEL OF ROTATING BOSE-EINSTEIN CONDENSATES 

ALEXEY CHESKIDOV, DANIEL MARAHRENS, AND CHRISTOF SPARBER 


Abstract. We study the long time behavior of solutions to a nonlinear par¬ 
tial differential equation arising in the description of trapped rotating Bose- 
Einstein condensates. The equation can be seen as a hybrid between the well- 
known nonlinear Schrodinger/Gross-Pitaevskii equation and the Ginzburg- 
Landau equation. We prove existence and uniqueness of global in-time so¬ 
lutions in the physical energy space and establish the existence of a global 
attractor within the associated dynamics. We also obtain basic structural 
properties of the attractor and an estimate on its Hausdorff and fractal di¬ 
mensions. As a by-product, we establish heat-kernel estimates on the linear 
part of the equation. 


1. Introduction 

1.1. Physical motivation. The study of quantized vortex dynamics in Bose- 
Einstein condensates (BECs) is a topic of intense experimental and theoretical 
investigations. A particular interesting situation is created when the BEG is stirred 
through an external rotating confinement potential. Indeed, if the rotation speed 
exceeds some critical value vortices and, more generally, vortex lattices are being 
created, see, e.g., mm for a broader introduction. 

From a mathematical point of view, rotating BECs can be described within the 
realm of a mean-field model, the so-called Gross-Pitaevskii equation m- In the 
following, we shall assume, without loss of generality, that the system rotates around 
the z-axis with a given speed 17 G R. Placing ourselves in the associated rotating 
reference frame, the corresponding mathematical model is a nonlinear Schrodinger 
equation (NLS) given by 

(1.1) idttp = —-A'0 -I- + V{x)'ip — 

Here, t G R, x G with d = 3, or d = 2, respectively. The latter corresponds to 
the assumption of homogeneity of the BEC along the z-axis (see, e.g., [71 H?], for 
a rigorous scaling limit from three to effective two-dimensional models for BEC). 
The parameter A ^ 0 describes the strength of the inter-particle interaction, which 
in this work is assumed to be repulsive. The potential V describes the magnetic 
trap and is usually taken in the form of a harmonic oscillator, i.e. 

(1.2) V(x) = w G R. 

Here, and in the following, we choose V to be rotationally symmetric for simplicity. 
All our results can be easily generalized to the case of an anisotropic harmonic 
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oscillator. Finally, ^ILip describes the rotation around the 2 ;-axis, where 

(1.3) L'll) := -i{xidx2i’ - X2dx^'ii)) = -ix-^ ■ Vi/;, 

denotes the corresponding quantum mechanical rotation operator. 

Most rigorous mathematical results on vortex creation are based on standing 
wave solutions of (ini,i.e. solutions of the form tp{t, x) = (p{x)e £ K, which 

leads to the following nonlinear elliptic equation 

(1.4) — 2 '^'^ + V{x)(p — VlLip — ^ip = 0. 

Equation (dll) can be interpreted as the Euler-Lagrange equation of the associated 
Gross-Pitaevskii energy functional El ED]: 

(1.5) Egp{p) ■■= dx, 

One possible way of constructing solutions to (11.41) is thus to minimize (11.51) under 
the constraint II 95 III 2 = M, where M > 0 denotes a given mass. This consequently 
yields a chemical potential p, = p{M) ^ 0 playing the role of a Lagrange multi¬ 
plier. In order to do so, one requires w > |0| which ensures that Eqp is bounded 
from below. Physically speaking, this condition means that the conhnement po¬ 
tential V{x) is stronger than the rotational forces, ensuring that the BEG stays 
trapped. Within this framework, it was proved in [30] that the hereby obtained 
physical ground states, i.e. energy minimizing solutions of dm), undergo a symme¬ 
try breaking (of the rotational symmetry) for sufficiently strong and/or A ^ 0. 
The latter is interpreted as the onset of vortex-lattice creation. 

On the other hand, it is often argued in the physics literature that a small amount 
of dissipation must be present for the experimental realization of stable vortex 
lattices, cf. dZlESlED]- In order to describe such dissipative effects, not present in 
the original Gross-Pitaevskii equation (du, the following phenomenological model 
has been proposed in [34] and subsequently been studied in, e.g., [TOl IT^ [25l [26] : 

(1.6) {ip — ^)dt'p = —-|- A|i/;pi/; -I- V{x)ip — ilLip — pip. 

Here /3 G M and 7 > 0 are physical parameters whose ratio describes the strength of 
the dissipation. (In |21) the authors use formal arguments based on quantum kinetic 
theory to obtain ^ « 0.03.) Note that any time-independent solution ip = <p{x) of 
dH]) solves the stationary NLS (II.4|) . In contrast to dn]), equation m is no 
longer Hamiltonian and only makes sense for t G R_|_. 


1.2. Mathematical setting and main result. This work is devoted to a rigorous 
mathematical analysis of dH])- In particular, we shall be interested in the long time 
behavior of its solutions as t ^ -l-oo. To this end, it is convenient to re-scale time 
such that /3^ -I- 7 ^ = 1. Then we can write 


iP-"f = -e'-^, 


for some id G 


( 


TT TTN 

2 ’ 27 ’ 


Note that by doing so, the real part of has the same (positive) sign as 7 > 0. We 
shall thus be concerned with the following initial value problem for (t, x) G M+ x 
and d = 2 ,3: 


(1.7) - e^^dtip =-^AipX\ip\'^'^ip-{-V{x)ip-nLip - pip, ipn^o = ipo{x). 


where ipo will be chosen in some appropriate function space (see below), and cr > 0 
a generalized nonlinearity. Formally, the usual Gross-Pitaevskii equation dm is 
obtained from dm in the limit id —±:|. On the other hand, if -d = 0 the 
Hamiltonian character of the model is completely lost and dm instead resembles 
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a nonlinear parabolic equation of complex Ginzburg-Landau (GL) type, cf. [5] for 
a review on this type of models. 

Equation (ini) can thus be seen as a hybrid between the Gross-Pitaevskii/Non- 
linear Schrodinger equation and the complex Ginzburg-Landau equation. Both kind 
of models have been extensively studied in the mathematical literature: For local 
and global well-posedness results on NFS, with or without quadratic potentials E, 
we refer to [HI [HIE]. Allowing for the inclusion of a rotation term, the initial value 
problem for dm has been analyzed in |H] . Similarly, well-posedness results for the 
complex GL equation in various spaces can be found in [191 US [2l] ■ The existence 
and basic properties of a global attractor for solutions to GL (on bounded domains 
U C are studied in m and [28] . Moreover, the so-called inviscid limit which 
links solutions of GL to solutions of NLS has been established in [SSj. However, 
none of the aforementioned results directly apply to the model (ini), which involves 
an unbounded (quadratic) potential V and a rotation term, neither of which have 
been included in the studies on GL cited above. One should also note that the GL 
equation in its most general form allows for different complex pre-factors in front of 
the Laplacian and the nonlinearity. In our case those pre-factors coincide, allowing 
for a closer connection to NLS. Very recently, a similar type of such restricted GL 
models with A < 0 (and without potential and rotation terms) has been studied in 
[T2l[T3] as an “intermediate step” between the NLS and the nonlinear heat equation. 
Finally, we also mention that equation with /3 = 0 is used to numerically obtain 
the Gross-Pitaevskii ground states, cf. [HidS]. 

As announced before, we shall mainly be interested in the long time behavior of 
solutions to ([13. In view of this the main result of our paper can be stated in the 
following form: 

Theorem 1.1. Let d G {2, 3}, oj > |H|, d G (—f, f), A ^ 0, and 0 < a < 2{d-2) ■ 
Then for any 

i/:oGE:={/GLri(M^) : \x\f G L^{R‘^)} 
there exits a unique strong solution ip G (^([O, oo), F) to (II. 7|) . The associated mass 
and energy thereby satisfy the identities (lOl) and 63 below. If, in addition, 
A > 0, the evolutionary system 63 possesses a global attractor A C S, i.e., A is 
is invariant under the time-evolution associated to ([13 and such that 

inf \\i}{t) - (^iiL2(Rd) 0. 

More precisely, 

A = {V'o : "00 = 0(0) for some 0 G C((—oo, oo); E) solution to 63 } 

is a connected compact set in L^(R^) and uniformly attracts bounded sets in L^(R^). 
Furthermore, for cr ^ A has finite Hausdorff and fractal dimensions which de¬ 
pend on the given parameters as described in Provosition 1 6.1 (A Finally, if p < ^ 
it holds A = {0}. 

Here, E is the physical energy space ensuring that AGp(0(t)) is finite. The 
assumption on ct > 0 is thereby slightly more restrictive than the one for the usual 
iL^-subcritical nonlinearities (see Remark |3.2I belowl. Note however, that we may 
always take cr = 1 in the above theorem which corresponds to the usual cubic 
nonlinearity. In addition, the condition lu > |H| ensures that the confinement is 
stronger than the rotation, and thus, the system remains trapped for all times t ^ 0 . 

As we shall see, neither the mass nor the (total) energy are conserved quantities 
of the time-evolution, but for A > 0, there are absorbing balls for 0 in both the 
mass and the energy space, see Section[5]for a precise definition. The existence of a 
global attractor A therefore requires the presence of the nonlinearity and, of course, 
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the presence of the confining potential V. Clearly, all stationary solutions ip € T, 
of (11.411 are members of A. However, since for fj, sufficiently large there are always 
at least two such solutions (namely, zero and the nontrivial energy minimizer) and 
since A is connected, it is unclear what the precise long-time behavior of dni) is. 
Indeed, in the case of the GL equation for superconducting materials it is known 
[55] that the global attractor contains not only all possible steady state solutions, 
but also the heteroclinic orbits joining these steady states, and we consequently 
expect a similar behavior to also hold also in our model. 

Except in the case /x < ^, the precise dependence of the dimension of A on 
the given physical parameters is not known. In Section [6.2l we shall prove that the 
Hausdorff dimension dimH(^) ^ m, where m depends in a rather complicated way 
on all the involved parameters. It is interesting, however, to check that m —>■ -boo, 
as |0| —>■ Lo. In other words, the influence of the rotation term potentially increases 
the dimension of the attractor. This is consistent with numerical and physical 
experiments on the creation of vortex lattices in rotating BEC. For a recent (non¬ 
rigorous) study which employs numerical simulations and asymptotic analysis to 
investigate the corresponding pattern formation mechanism, we refer to [lOj . In 
fact, one easily observes that in the linear case (A = 0) the dynamics admits expo¬ 
nentially growing modes, cf. Section [2T] below for more details. It is argued in [lO] 
that this type of instability mechanism is responsible for the nucleation of a large 
number of vortices at the periphery of the atomic cloud, as can be seen in physical 
experiments. 

The proof of Theorem 1 1.1 1 will be done in several steps: First, we shall establish 
local (in-time) well-posedness of (11.711 in Section |5| below. Then, we will show how 
to extend this result to global in-time solutions in Section [T] where we also prove 
that for ^ ^ solutions decay to zero as t -boo. The main technical step for 

the existence of an attractor is then to prove certain uniform bounds on the total 
mass and energy as done in Section 0 This will allow us to conclude the existence 
of an absorbing ball and of a global attractor in Section |6l where we shall also 
prove the announced estimates on the dimension. Finally, we collect some basic 
computations regarding the kernel of the linear semigroup in the appendix. 


2. Mathematical preliminaries 


In this section we shall collect several preliminary results to be used later on. 

2.1. Spectral properties of the linear Hamiltonian. In the following, we de¬ 
note by 



( 2 . 1 ) 


the linear Hamiltonian operator, with V{x) given in (11.21) . Note that in the case 
without rotation, i.e. = 0, the operator 



(-A+w^lip) , 


( 2 . 2 ) 


is nothing but the (isotropic) quantum mechanical harmonic oscillator in, respec¬ 
tively, d = 2, or 3 spatial dimensions. The spectral properties of Hq are well known 

[ini[55]: 

Lemma 2.1. Hq is essentially self-adjoint on C“(R‘^) C L^(R'^) with compact 
resolvent. The spectrum of Hg is given by a{Ho) = {Fo,„}neN; where 
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In addition, the eigenvalue -Eq,™ is ^ —fold degenerate. 

In particular, Eo,„ ^ Eo,i = ^ > 0, for all n G N. The associated eigenfunctions 
form a complete orthonormal basis of In d = 2, they are explicitly given 

by PO] : 

= fnAxi)fn2{x2), rij € N, 

where ni + n2 = n and the fn^ G iS(IR) are the eigenfunctions of the one-dimensional 
harmonic oscillator, i.e., an appropriately normalized Gaussians times a Hermite 
polynomial of order nj — 1. An analogous formula holds in d = 3 dimensions. 

In the case with 17 ^ 0, we first note that the commutator [Hq, L] = 0, due to the 
rotational symmetry of the potential V. This implies that Hq and L have a common 
orthonormal basis of eigenfunctions {XrajneNo: which can be obtained by taking 
appropriate linear combinations of the eigenvalues of Hq, see [20) . An important 
assumption throughout this work, will be that w > |r7|, ensuring confinement of 
the BEG. In mathematical terms, this condition implies that the rotational term 
can be seen as a perturbation of the positive definite operator Hq, such that Hq is 
still positive definite. In other words, we have that 

( 2 . 3 ) Hq\yi = Eqji)(ji, 

where the new eigenvalues S R (indexed in increasing order) are related to 

the unperturbed Eo,„ via 

{Eo^„, n G N} = {Eo,i + mlt, —I+l ^ m ^ i — 1, for£ G N}. 

In particular, under the assumption that w > 17, we still have: En^„ ^ for all 
n G N. Thus, the ground state energy eigenvalue stays the same with and without 
rotation. 

With these spectral data at hand, we can now look at the linear time-evolution 
(A = 0) associated to (11.611 . i.e. 

(2.4) {ij3 - "f)dt'tp = Hn'ijj - gLif- 

Using the fact that {xnjneN comprises an orthonormal basis, we can decompose 
the solution to this equation via 

(2.5) tpA, 2:) = ^ Cn{t)Xn{x), 

nSN 

where {cn(7)}„gN £ i-e. E |Cn(^)P < +00. In view of (12.3L (12.411 we find 
Cn{t) = c„(0) exp{-{iP + 'y){EQ,n - f-A), 

In particular, the normalization ,5^ -|- 7^ = 1 yields 

00 00 

Unimh ^ E = E 

n—1 n—1 

where we identify 7 = cosd. For d G (—f, f) the right hand side exponentially 
decays to zero as 7 —>■ -l-oo, provided /i < En,™, for all n G N. This is equivalent 
to saying that /i < Eq^i. On the other hand, if c„(0) = 0, then c„(7) = 0 for all 
7 > 0. Hence, given a // > Eq 1 the solution is exponentially decaying as long as 
the initial data is such that Cn(0) = 0 for all n G N for which En,ra < d- Otherwise, 
we have, in general, exponential growth of the E^-norm of '0(7). 

Remark 2.2. In the case where we choose /i = En,m for some fixed m G No, 
we see that the \cm{t)\^ = |cm(0)p is a conserved quantity of the linear time 
evolution. All higher modes exponentially decay towards zero, whereas all lower 
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modes will exponentially increase. We consequently expect linear instability of 
stationary states of the nonlinear system. 

2.2. Dispersive properties of the linear semi-group. In order to set up a well- 
posedness result for the nonlinear equation (HH), we need to study the regularizing 
properties of the linear semigroup associated to iJo, i.e. 

S'n(t) := exp t G R+, 

As usual we identify Sq {t) with its associated integral kernel via 

SQ{t)f{x)= [ Sn{t,x,y)f{y)dy, / G 

The following lemma states some basic properties of S'n(0 fo be used later on. 
Lemma 2.3. Let d G (—f, f) and t > 0. Then 

( 2 . 6 ) Sa(t.^.v) = 

where the pre-factor in front of the exponent is understood in terms of the principal 
value of the complex logarithm, and the phase function F is given by 

X, y) = cosh(e“*’’a;t) - cosh(e“*’’m)(a;iyi + X 2 y 2 ) 

+ i sinh(e”®’’f2t)(a;2yi - *12/2)^ ■ 

Moreover, for uj > |n|, there exists 6 > 0 such that 

(2.7) \\Sn{t)f\\L^^Cti^^-^.^\\f\\L, 

and 

(2.8) ||V5a(t)/||L^ + \\xSn{t)f\\L^, < C 

for all 1 ^ q ^ r ^ 00 and all 0 < t < 6, where the constants C and S only depend 
on d, Lo, and LI. 

The proof of Lemma 12.31 is a lengthy but straightforward calculation given in 
the Appendix. It is based on the well-known Mehler formula, cf. [3, and a time- 
dependent change of coordinates introduced in [4]. 

Remark 2.4. The decay estimates stated above are the same as for the heat 
equation. Indeed, Sq{T) may be viewed as an analytic perturbation of the classical 
heat semigroup. In the case without potential and without rotation, i.e. id = w = 0, 
similar estimates have been derived in, e.g., |12) . 

3. Local well-posedness 

In this section we set up a local well-posedness result for the initial value problem 
(ini). In order to do so, we use Duhamel’s formula to rewrite (EB as 

(3.1) '0(t) = 5a(t)'0o - / Sni^t - T){X\tj}{T)f'^ - p)'tj}{T) dr, 

Jo 

for all t G [0,r]. Here, and in the following, we denote ^^{t) = '0(t, •). We shall 
work in the physical energy space given by 

E = {f : \x\f G L^{R‘^)} 

and equipped with the norm 

ll/ll|:= ll/lli^ + llV/lli^ + llx/lli.. 




GLOBAL ATTRACTOR IN A MODEL FOR ROTATING BEG 


7 


The estimates on the semi-group Sn(t) stated in Lemma 12.31 allow us to infer the 
following result (which is similar to those in [T21 [23] ) • 

Proposition 3.1. Let A,/! G K, i? G j f); ^ ^ ^ {2)3}. 

(i) Let p > max((Td, 2(7 -|- 1) and tpo G LP{W^). Then there exists a time 

T > 0 and a unique solution ip G (^([O, T]; to (11.71) . depending 

continuously on the initial data. 

(ii) If, in addition, 0 < a < 2 (d- 2 ) V'o £ then there exists aT* > 0 such 
that the solution from (i) satisfies 

iAGC([0,n;E). 

Moreover, the solution is maximal in the sense that either T* = -foo, or 
the following blow-up alternative holds: 

lim ||i/;(t)||E = oo. 

+ _ i,T'* 


Proof. The proof is based on a fixed point argument using Duhamel’s formula and 
the properties of the semigroup 5'n(<). To this end, we first note that the term pijj 
is of no importance here, as it can always be added in a subsequent step (in fact, 
we could have included it in the kernel of Sn{t)). Hence let us assume that p = 0 
for notational convenience. 

To prove (i), we will show that the mapping 


'll; := Sn{t)iljo - e [ Sn{t - dr 

Jo 

is a contraction in the space 

XT:={fj€C{[0,T];LP{R‘^)) : ML^io.T-L.) ^ 

for small enough T > 0. To do so, we can use the kernel estimate ILZI) with the 
following choice of parameters: 


r-p^2a+l. 

when 

II 

P 

p > max(crd,2cr-f 1), 9 = 2cr + 1 ’ 

when 

d = 3. 


Note that any such a choice of p implies that da < p. One can also see that 
l^g^r^ooin both cases. This yields 

W^WiOhp < UoWlp + a / \\Sn(t - t){\iP{t)\^'^iP{t)) dr 

Jo 

^ WMlp + c fit - r)-"'^/^||i/>(r)||2‘}+i dr 
Jo 

^ \mLp + I dr 

Since cr < |, the remaining integral is finite and hence, 

mm)\\LP ^ + 
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where 1 — ^ > 0. Thus, for T > 0 sufficiently small, we conclude that ^ indeed 
maps Xt onto itself. Likewise it holds that for two solutions '0 and 0 

||s( 0 )(t)-s( 0 )(t)iu. 

< A / ||S'n(t-r)(|0(r)p'^0(r) - |0(r)p'^0(T))||^p dr 
^0 

0 C r(i-r)-^‘^/^(||0(T)||r. + Il0(r)|ir.)||0(r) -0(r)|U. dr 

0 CT P ||0||£So('o^T;L!>)II'*^ ~ ''/’IU“(0,T;LP)i 
which shows that S is a contraction for T > 0 sufficiently small. 

To prove (ii), we first note that by Sobolev imbedding E ^ LP(R^), for p < 
p* = when d = 3 and p < oo when d = 2, respectively. Thus E n LP(W^) = E 
for p < p*. We now want to show that for 0 < cr < 2 (d- 2 ) ’ ^ norm of the 

solution is controlled by an appropriately chosen Lp norm satisfying p < p* and the 
conditions in part (i). 

The first step to do so, relies on appropriate expressions for the commutators 
[V, S'n(t)] and [a;, S'n(t)]. At least formally, it holds that 

-e*’’at[V,5a(t)] = [V,ida^nW] = Hn[X,Sn{t)] + [V, ida]5n(t), 


and one easily computes 

[V, Hn] = vy + 


in view of Hence, Duhamel’s formula and the fact that [V, 5*0(0)] = 0 imply 

[V,5'n(t)] =-e”®’’ f Sn{t - t)(uj'^x + irty-'-)SQ{T) dr. 

Jo 


Using the fact that 
we likewise obtain 


[x, Hq] = V — iUa;'* 


[x, 5'n(t)] =-e f - T){y - ii}x-'-)Sn{T) dr. 

Jo 

Straightforward calculations then yield 

V0(t) = 5n(i)V0o - e-®’’ [ Sa{t - r) (AV(|0|2-0) + + zL!V^)0) (r) dr 

do 

as well as 

xip{t) = S'n(t)a;0o - e“®’’ f Sn{t - r) (Aa::|0p‘^0 + (V - *Ha:-‘-)0) (r) dr. 

do 

We consequently expect that the combination of 0,x0, V0 will form a closed set 
of estimates (a fact already observed in [4]). 

It follows that the E norm of 0 is controlled by its LP norm. For instance, choose 
r = 2 and q such that 

1 2cr 1 . ,_, 

m ([2Jl). 

Note that the condition g 0 1 is equivalent to p Aa and thus we require the 
existence of a p such that 


2d 


max(4(T, ad, 2cr + 1) < p < 


d-2' 
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For d £ {2,3} one can readily check that 0 < cr < 2 (d- 2 ) ensures the existence of 
such a p. It follows that 

f S'a(t-r)V(|V'P‘"V’)(r) dr 

Jo ^ 

pt 

< C'll^lli~( 0 .T;LP)l|VV’IU“( 0 ,T;L 2 ) / dr. 

Jo 

On the other hand, the linear terms can be estimated with r = g = 2 in (lO) , to 
obtain 

IIVV^IU“(0,T;L2) < ||V^o||l2 + P IIV’IIl~(0,T;Lp)II^V^IIl“(0,T;L2) 

+ C 2 T (||a;'l/'||L=x=(o,T;L2) + ||Vi/)||loo(o,T;L2)) 

< IIvv^oIIl^ + 

for r > 0 sufficiently small (depending on the size of ||'i/'||l”Lp)- Similar arguments 
for ip and xtp imply 

IIV'IIl“(0,T;S) ^ IIV'oIIe + CT^ p IIV'lli~(0,T;LP)IIV'llL“(0,T;E)- 

Choosing T > 0 even smaller, if necessary, the second term on the right hand side 
can be absorbed on the left hand side and we are done. As before, this inequality 
also applies to the differences of two solutions ip, ip, which yields the continuity of 
Ip in E. 

We denote by T* > 0 the maximal time of existence in E. This is always less 
than or equal to T > 0, the maximal time of existence in Lp{W^). To prove the 
blow-up alternative, assume by contradiction that T* < oo, and \\ip{t, •)||e remains 
bounded for t £ [0,T*]. Then, by Sobolev imbedding ||i/'(t, OIIlp also remains 
bounded and thus, we can restart the local existence argument in E leading to a 
contradiction. □ 

Remark 3.2. Unfortunately, our method of proof does not yield existence of solu¬ 
tions for the full iJ^-subcritical regime, i.e., a < We expect that this is only a 
technical issue that can be overcome using a different approach (for example, by us¬ 
ing ideas from [^, or by generalizing the space-time estimates of [B] to Sq). Note, 
however, that our slightly more restrictive condition a < 2 [d- 2 ) allows to take 
CT = 1 in d = 3. Hence, the physically most relevant case of a cubic nonlinearity is 
covered. 

4. Global existence and asymptotic vanishing of solutions 

In this section, we shall first prove the global existence of solutions in the energy 
space before showing that for any choice oi p < Eq^i, the solutions asymptotically 
vanish as f —>■ -l-oo. 

4.1. Global existence. In order to prove global well-posedness of dni), we will 
need to collect some useful a-priori estimates. To this end, we denote for ip £ E 
the total mass by 

(4.1) M{iP) := WiPWl., 
and the total energy by 

(4.2) E{iP) := [ (hvipf + V{x)\iPf + - nPpLip] dx. 

The latter is nothing but the sum of the kinetic, potential, nonlinear potential, 
and rotational energy. Clearly, for ip £ T,, Sobolev’s imbedding implies that all the 
terms in E{ip) are finite, provided a < (and hence also for our range of a). 
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For simplicity of notation, we will write E{t) = •)) and likewise for M{t), 

whenever we compute the mass and energy of the time-dependent solution tp(t,x) 
to (dH). In addition, the free energy is given by 

(4.3) Fiif) := Eif;) - 

In the case of the usual Gross-Pitaevskii equation, i.e. i? = ij, one finds, that 
both M{t) = M(0) and E{t) = E{0) are conserved in time [3]. In our dissipative 
model this is no longer the case. Instead we have the following result, which can 
be seen as an extension of some well-known identities proved for the classical GL 
equation, cf. [mu El [35]. 

Lemma 4.1. Let a < 2 (^- 2 ) ^ ^^([0, P]; S) be a solution to (II.7|) . Then the 

following identities hold: 

(4.4) M{t) + 2cosi?£ (^E{s) + - /iM(s)^ ds = M(0), 

and 

(4.5) F(t)-|-2cosd f f \dt4’is,x)\'^ dx ds = F{0). 

Jo Js.-^ 

In particular, for D G (—f, f), the free energy F['if) is a non-increasing functional 
along solutions of (ini). 

Proof. In a first step, let us assume sufficient regularity (and spatial decay) of ip, 
such that all the following calculations are justified. Then, as in the case of the 
usual NLS, identity (14.41) is obtained by multiplying (11.71) by ip, integrating with 
respect to x G and taking the real part of the resulting expression (see, e.g., 
mn]). This yields 

(4.6) |M(f) = -2cosd (^E{t) + - pM{t)^ 

which directly implies ( 031 ) after an integration in time. Similarly, after multiplying 
( 031 ) by dtip, integrating with respect to x, and taking the real part, we obtain 

d f 

(4.7) -—[E{t) — fj,M{t)) = —200511 / \dtipit,x)\^ dx, 

dt J^d 

which yields ( 031 ) after integration w.r.t. time. 

The second step then consists of a classical density argument (cf. 02]), which, 
together with the fact that ip(t) depends continuously on the initial data ipo G E, 
allows us to extend (|4.4I) and (14.51) to the case of general solutions ip G (^([O, T; E). 
Finally, we note that for 1 } G (—f, f) we have cosd > 0, and thus (1331) directly 
implies that F(t) ^ F(0), for all t ^ 0. □ 

Having in mind that ip G C'([0,r], E) the assumption on a implies (via Sobolev 
imbedding) that the integrand appearing in identity (14.41) is a continuous function 
of time. The fundamental theorem of calculus therefore allows us to differentiate 
(14.41) w.r.t. t and consequently use the differential inequality (14.61) . However, the 
same is not true for (1331) . i.e., we cannot use (1331) . since at this point we do not 
know wether dtip G (^([O, T; L^(IR‘^)) holds true. This fact will play a role in some 
of the proofs given below. 

Another preliminary result, to be used several times in the following, is the 
fact that under our assumptions on the parameters w, H, A, tr, the energy is indeed 
non-negative. 
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Lemma 4.2. Let ui > |f2|, A ^ 0, and a < Then for any u G T, there exists a 

constant c = c(uj, LI, X, a) > 0, such that such that 

\\S/u\\l. + \\xu\\l. + \\u\\%t^.^cE{u). 

Proof. Since A ^ 0, the only possibly negative term within E(u) is given by the 
rotational energy. However, since Lif jitp' =: e < 1, Young’s inequality applied 
to (O yields the pointwise interpolation estimate 
2 

|OuLu| ^ dx + ^ V{x)\u\^ + ^|Vup. 

We therefore can bound the energy from below via 

0 ^ ^llVulli. + ^ E{u). 

Analogously, we have 

0 ^ ^ E{u). 

Combining these two estimates then yields the desired result with a constant 

4 

min{l-e,|^}' 

Note that c —)■ +oo as |H| -g- oj. □ 

The mass/energy-relations stated in Lemma iTlI can now be used to infer global 
existence of solutions in the case of defocusing case A > 0. 

Proposition 4.3. Let w > |H|, i? G A ^ 0, and a < 2 {d- 2 ) ■ Then, for 

any i/iq G S there exists a unique global-in-time solution tp G C([0,oo); E) to (11.711 . 

Proof. In view of the blow-up alternative stated in Proposition 13.11 all we need to 
show is that the E-norm remains bounded for all t ^ 0. Lemma [4.21 implies that 
this is the case, as soon as we we can show that both M{t) and E{t) are bounded. 
In order to do so, we first consider the case /i < 0 and recall that cosi? > 0 for 
I? G (—§, §). In this case identity (14.51) implies 

E{t) -\- \^\M{t) < F(0) < -1-00, 

and since both E(t) and M(t) are non-negative, we directly infer the required bound 
on the mass and energy. 

On the other hand, for /i ^ 0, identity (lOl) yields (since A > 0) 

M(t) < M(0)-I-2^cosd f M{s)ds, 

Jo 

and hence, Grownwall’s lemma implies 

(4.8) M{t) ^ M(0) (1 -h 2/rtcosi?e2'"*™"’’) . 

Using this estimate in identity (H31) we obtain 

E{t) < F(0) -h /iM(t) < £;(0) -f 2/i2tcosi?M(0)e2'"*™"’’. 

The right hand side is finite, for all t ^ 0 and thus, the assertion is proved. □ 

Remark 4.4. The global in-time strong solutions constructed above are of the same 
type as the corresponding solutions for NLS with quadratic potentials, cf. BE- 
It is certainly possible to, alternatively, construct global weak solutions to (ini) 
as has been done for the usual GL model in, e.g., [TH [23]. But since we consider 
the equation (II3 as a toy model describing possible relaxation phenomena in the 
mean-field dynamics of BEC, we have decided to remain as close as possible to the 
corresponding NLS theory. In particular, we do not make any use of the strong 
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smoothing property of the linear (heat type) semigroup Sn{t) for G We 

hnally note that our set-up makes it possible to directly generalize the inviscid limit 
results of [35] to our model. 

4.2. Asymptotically vanishing solutions. The discussion in Section |2T] shows 
that solutions to the linear time evolution A = 0 asymptotically vanish, provided 
< Eq, i.e., the lowest (positive) energy eigenvalue of Hq. We shall prove that the 
same is true for in the nonlinear case A > 0. 

Proposition 4.5. Let d G (—f,f); A > 0, a; > |f2|, and ijj G (^([OjOo),^) be a 
solution of (IL3 with fi < i?n,i = Then 

lim ||V'(t)||L2 = 0, 

i—>-+oo 

exponentially fast. 

Proof. For solutions if G (^^([O, oo), E) we are allowed to use the differential inequal¬ 
ity (ITO . which together with the fact that A ^ 0 implies 

^M(t) ^ —2 cos'd {E{t) — pM{t)). 

Decomposing ip(t, x) in the form (12.511 . and dropping the nonlinear term oc A within 
E(t), then allows us to rewrite this inequality as 

d ^ 

—M(t) < -2cos (Fin,„ - M)|cn(i)P < -2cosi?(£ln,o - d)M{t), 

n—0 

since Fln.n — M ^ Fln,o — h- > 0, and M{t) = \^n{t)\‘^. The inequality above 

can thus be rewritten as 

^ ^^+2tcosi)iEa.o-u)M{t)'^ 0 , 

which after an integration in time implies 

since 'd G (—f, f )■ □ 

At this point, it is unclear if the decay rate given above is indeed sharp. 

Remark 4.6. In the case where // < 0, one does not need to use the decomposition 
of if via the spectral subspaces of Tin , at the expense of a slightly worse decay rate. 
Indeed, for /r < 0, the inequality (14.61) directly yields 

^M(t) ^ —2\p\cos'dM{t), 

and thus 

M{t) s: M(0)e-2*l'"lVt ^ 0. 

Note that for /r < 0 there are no nontrivial steady states (p{x) ^ 0, satisfying (II. 4|) . 
This can be seen by multiplying equation CH) with (/3, integrating in x G R'*, and 
recalling the restriction w > D ^ 0, which implies that p, has to be non-negative. 

5. Bounds on the mass and energy 

In this section we shall prove the existence of absorbing balls in both L^(R‘^) 
and E for solutions to (113. In view of the discussion on the linear model, cf. 
Section \TJ\ this might seem surprising, given that for general ^ > 0 we can expect 
exponentially growing modes. However, we shall see that for A > 0, the nonlinearity, 
combined with the confining potential, mixes the dynamics in a way that makes 
it possible to infer a uniform bound on the mass and energy, and hence on the 
E-norm of the solution. To this end, the following lemma is the key technical step. 
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Lemma 5.1. Let A > 0, w > |f2| and 0 < ct < 2 (d- 2 ) • Then there exists a constant 
C = C{u}, O, A, a) > 0, such that 

^ , with d= . 

Proof. The proof of this result relies on the following localization property: For all 
d ^ 1 and all p ^ 2 and any / G 

(5-1) ll/llL2(Rd) ^ 2||a;/|| j;^2(Rd) 

with 

g ^ ~ p) ^ d{p - 2) 

l + d(i-i) 2p + d(p-2)- 

In order to show this, let Bn denote the ball around the origin of radius r > 0. We 
rewrite 

ll/llL2(R<i) = ||/||l2(b^) + ||/||L2(Rd\B^) ^ r‘^^^ ’’^||/||lp(Bp) + -||ai/||L2(Rd\B^) 

^ ||/||^p(Rd) + -||a;/||i2(Rd). 

The right-hand side is minimal if both summands are of the same order of magni¬ 
tude, i.e. 

^l+d(i-i) ^ WxfWL^jRd) 

ll/llLP(R<i) 

With this choice of r, the estimate (EU follows and a density argument allows to 
extend it to any / G S. Specifying p = 2a + 2, consequently yields 

(5.2) IIV^IIi2 ^ 2(" / \x\'^\iIj{x)\'^ dx'\ ( [ \^{x)\^^+^ dx) 

\jR‘i /V Js.'i J 

where 9 = 2 a+ 2 +da ■ view of Lemma 14.21 both factors on the right hand side of 
(EH) are bounded by the energy. More precisely, 

M{if) ^ 2{cE{i})Y+^ = CE{if)^, 

a9 + l 

where C = 2c ‘'+1 and c = c{uj, ft, X,a) >0 is the constant from Lemma 14.21 □ 

Remark 5.2. Note that in order to infer this bound one needs the presence of both 
the confinement and the nonlinearity, since the proof requires cr > 0, A > 0 and 
w > 0. Moreover, one checks that C —-boo, as |fl| —>■ oj. 

With this result in hand, we can deduce global bounds on M{t) and E{t) along 
solutions of E3- 

Proposition 5.3. Let -0 G (^([O, oo),S) he a solution to p.7l) with 1? G (—f,f)- 
Under the assumptions of Lemma \5.1l if additionally fi> 0, there exists a constant 
K = K(uj, 17, a, A, p) > 0, independent of time, such that 

E{t) sf K+ e-*f^^°^^E{0), Vt^O. 

Proof. We first note that Lemma 15.II and the differential inequality (14.61) imply 

^M{t) ^ —2 cost) Eit) + C pE{t)~^. 
dt 

Now, for any d G (—ff) and 9 = by Young’s inequality, we obtain 
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where C > 0, depends on all the parameters involved, but not on time. Thus, we 
have 

^ — cos ’OElt) + ^lCC. 
at 

On the other hand, identity (1131) implies 

E{t) — E{to) < ^M{t) — 0<to<t, 

and hence 

E{t) — E{s) ^ f {—ficosd E{t) + fi^CC) dr, 0 < s < t, 

J S 


as well as 


E{t) — E{s) ^ j cos'd E{t) + fi^CC) dr, 0<t 

J S 


< S. 


Now, given any positive bump function x G — €,t + ej), such that x' ^ 0 

on (t — e,t) and x^ ^ 0 on + e), we multiply by x^(s) and integrate in s, to 
obtain 


/•t+e 

u 


[E{t) - £'(s)]x'(s) ds ^ 


nt-\-€. nt 


{—fj, cosd E{t) + fi'^CC)x'{s) dr ds 


-I 


t — e J s 

t pT 


{—fi cosd E{t) + fi^CC)x'{s) dsdr 


t — eJ t — e 

t-\-e f‘t-\-e 


{—fi cosd E{t) + fi^CC)x'{s) dsdr 


pt-\-e 


{—cos'd E{t) + ^^CC)x{t) dr. 


A similar computation gives the same inequality for a negative bump function 
function x € {{t —t,t +e)), such that x^ ^ 0 on (t — e, t) and x' ^ 0 on (t,t + e). 

Since an arbitrary test function can be written as a linear combination of positive 
and negative bump functions, we have 

- [ E{T)x'iT)dT f (-noos'dE{ t) + ^l‘^Co'] xi'r)dT, 

Jtn Jtn ^ ' 


for any x G t)). Here, we have also used the fact that x has compact 

support on (to,t). Choosing x(t) = ^(t){T) with (f> G we obtain 

-[ E{T){e^^^°^^(i){T))' dr [ [-ncosdE{T)+n^cd) 4 ){t) dr, 

Jtn Jtn ^ ' 


'to 

and thus 


— / E{T)e 


flT COS / 


0 '(r)dr< / n Hr) dr 


'to 


/ 

Jtc 

f dr. 

Jtn d cos'd 


Hence 


ii‘^r<n ii^EC 

^ (1 - ^ (1 - 

n cos'd n cos'd 

for almost all 0 < to < t. In summary, for almost all t ^ 0 we have 

E{t) < + A:(1 - 


K = 


ncc 

cos'd 


where 
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However, since '0 € ^([0, cx); E)) implies that E(t) is continuous in time, we conse¬ 
quently infer the inequality for all t ^ 0 . □ 

Remark 5.4. The proof above is slightly complicated due to the fact that we 
cannot use the energy identity (EH) in its differentiated form (I4T1) . see the discussion 
below the proof of Lemma |4.II If we ignore this problem for the moment, then we 
have 

^E{t) ^ ^ —fj,cos'dE{t) + 

Lc L lAj U 

which directly allows us to conclude the assertion proved above. 

In view of Lemma lS.li the bound on E(t) obtained above implies a similar bound 
on M[t). In particular, there is some constant pM > 0 and a function such 

that for all "0 G (^([O, oo); S) solutions to (II. 7L it holds 

Il0(i, OIIl^ ^ PM, Vt>tM(M(0)). 

Therefore 

{0 G L2(M^) : 11011^2 0 pm} C 

is an absorbing ball for trajectories t i—> 0(t, •). Similarly, we know, that there 
exists a pe 0 pm and a function tE(-), such that 

llV'Ci, Oils 0 PE, vt 0 tE(||0(O)||E)- 

In other words, 

(5.3) X:={0gE:||0||e<Pe} 

is an absorbing ball in E for trajectories t i—>■ 0(0 •). In our study of long time 
dynamics of dni), the set X will play the role of a phase space. 


6. The global attractor and its properties 


In the previous section we proved that solutions ip(t) exist globally in E, and, 
moreover, all such solutions remain within an absorbing ball X G T, for t > 0 
large enough. It is therefore natural to ask whether there exists an .4 C E that 
attracts all trajectories t i—>• ijj{t, •) G E. Unfortunately, classical theories of global 
attractors (see, e.g., [ 111121 ]) do not apply to our situation as they typically require 
asymptotic compactness, which is unknown in E. However, the trajectories might 
still converge to the global attractor A in some weaker metric, say L^. To prove 
this we revisit the rather general framework of evolutionary systems introduced in 
m and adapt it to our situation. 


6.1. Existence of a global attractor. First, recall that our phase space is the 
metric space (X, di 2 (•, •)) where X C E is given by (15.31) and 6.^2 ( 0 , 0 ) = ||0 — 0 || ^^2 . 
We note that X is dL 2 -compact. In addition, we also have the stronger E-metric 
dE(0,0) := 110 — 0||e on X, which satisHes: If dE(0n,0n) ^ 0 as n oo for some 
■071, 0n G X, then dj;,2(0„,0„) —J, 0 as n oo. Note that any E-compact set is 
L^-compact, and any L^-closed set is E-closed. 

Now, let C'([a, 6]; X,), where • = E or be the space of d,-continuous X-valued 
functions on [a, b] endowed with the metric 

f^C([a.b];A.)('0,0) := sup d. (0(t), 0(t)). 

[a,b] 


Also, let C{[a, oo); X,) be the space of d,-continuous X-valued functions on [a, oo) 
endowed with the metric 


'^C([a.oo);A.)(0, 0 ) 


E 

TeN 


1 sup{d,(0(t), 0(t)) : a < t < a + T} 

2^ 1 -I- sup{d,(0(t), 0(t)) : a < t < a + T} 
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In order to define a general evolutionary system, we introduce 
T := {/ : I = [T, oo) C M, or I = (—00, 00)}, 
and for each I C T, we denote the set of all X-valued functions on / by X{I). 

Definition 6.1. A map £ that associates to each / g T a subset £{I) C X{I) will 
be called an evolutionary system if the following conditions are satisfied: 

(i) £:([0,oo)) ^ 0. 

(ii) £{l + s) = {i/'(-) ; tj}{' + s) £ £{!)} for all s £ R. 

(iii) For all pairs I2 C h gT: {V'(')l-f2 ^ ''/'(•) £ C Sih)- 

(iv) £:((-oo,oo)) = {!/>(•) : '(/'(•)I [t.oo) £ £{[T,oo)) VT £ R}. 

In general, £{I) will be referred to as set of trajectories on the time interval /, and 
trajectories in £{{—00,00)) will be called complete. 

We now consider the specific evolutionary system induced by the family of tra¬ 
jectories of inn in X. More precisely, we set 

£{[T, 00)) := I '0 £ C!{[T, 00); X) a solution to (I1.7I1 . with i? £ ( — 

(6-1) d 1 

A,/r > 0, w > |fl|, and 0 < CT < 2 (^d- 2 ) \' 

Clearly, the properties (i)-(iv) above hold for the evolutionary system associated to 
(ini). In addition, due to Proposition l 4 . 3 l for any ifo G X there exists tp £ £{[T, 00)) 
with ■ 0 (r) = ipQ. Standard techniques then imply the following lemma: 

Lemma 6.2. Let {'ipn)neti be a sequence of functions, such that xpn G f([Ti,oo)) 
for all n £ N. Then for any T2 > Ti there exists a sub-sequence {ipuj ) jeN which 
converges in C([Ti, T2]; ^^2) to ip € f([Ti,oo)). 

Proof. Since X is compact in L^(R‘^), there exists a sequence {'tpnj)jGN such that 
'<pnj{Ti) —>• Ip for some ip £ L^(R'^). However, since lower-semicontinuity and the 
definition of X yield 

||0||s < liminf ||0„ (ri,-)||s ^ ps, 

j->oo 

we have that ip G X. In view of proposition 14.31 there exists ip £ f([Ti,oo)) with 
ip{Ti) = Ip. Continuous dependence on the initial data, then gives the desired 
result. □ 

Using this, we can prove one of the main structural properties of the set of 
trajectories induced by (II 3 : 

Proposition 6.3. £’([ 0 , 00)) is a compact set in (^([O, 00); Xi2). 

Proof. First note that f([ 0 , 00)) C (^([O, 00); Ari2). Now take any sequence £ 

f([0, 00)). Thanks to Lemma [521 there exists a subsequence, still denoted by ipn, 
that converges to some ip^ £ £'([ 0 , 00)) in (^([ 0 , 1 ]; ^^2) as n —00. Passing to 
a subsequence and dropping a subindex once more, we obtain that ipn —t ip'^ in 
C([ 0 , 2];Xl2) as n — >• 00 for some ip"^ £ £^([ 0 , 00)). Note that ip^{t) = ip'^{t) on [ 0 , 1 ]. 
Continuing and picking a diagonal sequence, we obtain a subsequence ipn^ of ipn 
that converges to some ip £ £"([0, 00)) in (^([O, 00); ^^2) as Uj 00. □ 

In order to proceed further, we denote, as usual, the set of all subsets of X by 
P{X). For every t ^ 0 , we can then define a map R{t) : P{X) P{X), by 

R{t)A := {ip{t) : ip{ 0 ) £ A, such that ip £ £([ 0 , 00))}, for any A C X. 

Note that the assumptions on £ imply that R{s) enjoys the following property: 

( 6 . 2 ) R{ts)A c R{t)R{s)A, AcX, t,s^ 0 . 
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Definition 6.4. A set A is called invariant under the dynamics, if R{t)A = A for 
all t ^ 0. 

We also recall the standard notion of and w-limit associated to an evolutionary 
system (see also [?T]1. 

Definition 6.5. The uj,-limit (• = E, of a set A C A is 

uj,{A) := f| U R{t)A. 

TJsO 

We also note that an equivalent definition of the w,-limit set is given by 
w,(A) = {tjj G X : there exist sequences oo and G R{tn)A, 

such that tpnitn) 4’ in the d,-metric}. 

Finally, we will give a precise definition of what we mean by an attractor. 

Definition 6.6. A set A C A is a d,-attracting set, if it uniformly attracts A in 
d,-metric, i.e. 

liminfd,(i?(t)A,((i) Q. 

ifi&A 

A set A C A is a d,-global attractor if A is a minimal d,-closed d,-attracting set. 
After these preparations, we are able to prove the main result of this section: 

Corollary 6.7. The evolutionary system (|6.1I) possesses a unique d]^ 2 -global at¬ 
tractor A = 072,2(A), which has the following structure 

A = {'00 : "00 = '0(0) for some if G £{{—oo, oo))} 

Furthermore, it holds: 

(1) For any e > 0 and T > 0, there exists a to € K, such that for any t* > tg, 
every trajectory ip G £"([0,00)) satisfies d2,2(0(t), 0(t)) < e, for all t G 
[t*,t* -\-T], where 0 G £1((—oo,oo)) is some complete trajectory, i.e., the 
uniform tracking property holds. 

(2) If the S global attractor exists, then it coincides with A. 

(3) A is connected in L^. 

(4) A is the maximal invariant set. 

Proof. Assertion (1) and (2) follow from the results proved in [IS]. To this end, 
one first shows that the o;2,2-limit of A is an attracting set, which by definition is 
closed and the minimal set satisfying these two properties. Then, using Proposition 
IQ and a diagonalization process, one can prove the structural properties of A, cf. 
[T51 Theorem 5.6]. The fact that A is connected then follows from Lemma 16.21 
and uniqueness: We argue by contradiction and hence assume that A is not T^- 
connected. Then there exist disjoint di2-open sets Ui,U 2 G A such that A C 
Ui U U 2 and An Ui, A n 1/2 are nonempty. Define 

X,={ifGX:uJL2{iP)GUj}, J = l,2. 

Since Ui and U 2 are disjoint, we also have that Ai, A2 are disjoint. Continuity 
of trajectories implies that Ai U A2 = A. Since A is di2-attracting, there exists 
T > 0 such that 

i?(t)A G C/i UC/2, \/t>T. 

By continuity of trajectories we have that for each ip G f ([0, 00)), either ip{t) G Ui 
for all t > T, or ip{t) G U 2 for all t > T. This implies that both Ai and A2 
are nonempty. Moreover, Lemma |6.2I implies that Ai and A2 are d2,2-open. This 
contradicts the fact that A is d2,2-connected. Finally we note that the structure of 
A, together with uniqueness of solutions, imply that A is an invariant set. Clearly, 
only complete trajectories are invariant, hence A is the maximal invariant set. □ 
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Remark 6 .8. In the case of the usual GL equation (posed on bounded domains 
D C K.'^) many more details concerning the global attractor are known, see, e.g., 
[28l [3n [33] . It is an interesting open problem to check which of these results can be 
extended to our situation and what the main structural differences between dm) 
and the usual GL equation are. 

6.2. Dimension of the attractor. We hereby follow the, by now, classical theory 
of estimating the Lyapunov numbers associated to f ([0, oo)) by studying the evolu¬ 
tion of an m-dimensional volume element of our phase space X, cf. [5n Chapter V] 
for a general introduction. Using this technique, the case of the usual GL equation 
on bounded domains D C K", with n = 1,2 is studied, e.g., in |3TJ Chapter VI, 
Section 7]. In our case, the same idea works, but requires several adaptions on a 
technical level. 

To this end, we first rewrite dni) as 

dtip = tjj^t=o = ' 00 - 

and, for any i/jq £ .4, consider the linearization around a given orbit 
i.e., 

(6.3) = V'|t=o = e- 

Here, ^ G X and G' denotes the Frechet derivative 

G'{'tp)cj) = Hq(I> - ficj) + X {tijcj))) , 

where Hq is the linear Hamiltonian (with rotation) defined in (12.11) . It is easy 
to see, that the linearized equation (16.31) admits a unique strong solution for any 
given ^ G X and tp G A. We now consider ,..., (pm(t) solutions to (16.31) , 
corresponding to initial data ■ Arm m £ N, and choose an L^-orthonormal 
basis - ■ ■ ,Xm{t) of 

Pm(t)X := span{(l)i{t),. ■ ■ 

where Pm denotes the corresponding orthogonal projection. Then, it is easy to see 
(cf. [31)1. that the evolution of the m-dimensional volume element in X is given by 

\(j)i{t) A • • • A <(<m(t)| = ICi A • • • A exp J ReTre“*’’G'(V'(s)) o PmA) ds 

In order to proceed, we first note that: 

Lemma 6.9. Let Hq he given by (12.21) . Then, for any orthonormal family {xj}jLi C 
L^(]R‘^) there exists a constant c = c{uj, d) > 0, such that 

m 

1=1 

Proof. Having in mind the form and the degeneracy of the eigenvalues stated 
Lemma [23 one checks that when counted with multiplicity Eo^m as 

m —1 oo. The desired result then follows directly from EH Chapter VI, Lemma 

2 . 1 ]. □ 

Using this, we can prove the following result for the dimension of A: 

Proposition 6.10. Consider the dynamical system (16.11) with c ^ and letm G N 
be defined by 

. X d/(d+i) 

m — 1 < ( — ^ 
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where 



K2 = CJfJ, 


1+d 



— d 


c"(A|^|)i+“ 



with c,c',c",a,a positive constants depending only on (jj,d,a, and 

rt 

II td! „\^l. I, 

IL2<t + 2 ' 


5 = lim sup sup I - 

t^oo tjjoeA \ t 


1 


||i?(s)i/>o||??f+2ds U ( iC- 


/o 


cr + 1 \ 


A 


Here, K is the constant from Proposition 15.^ 

Then, as t ^ +oo, the m-dimensional volume element in X is exponentially 
decaying. Moreover, the fractal (and hence Hausdorff) dimension of A is less than 
or equal to m. 


Proof. Having in mind the representation formula for the m-dimensional volume 
element as given above, we introduce 

gm := limsup sup f HeTr e~^^G'{ip{s)) o P^(s) ds\ . 

ll«dll,2=Sl V t Jo ) 

and quote the following result from [H Chapter III, Corollary 4.2]: If there are 
constants ki 2 ^ 0, such that 

fo ^-Kl/ + K2, Vj > 1 , 

then the fractal dimension of A enjoys the following bound: 

/ \ V® 

In order to obtain the required estimate on qj , we first note that 

m 

ReTre-*’’G'(i/>(t)) oP„(t) = ^Re(e-*’’G"(iA(t))x, (<),X,(t))i2. 

1=1 

Next, we recall that = ^ + i/d, with 7 > 0, and compute (suppressing all the 
t-dependence for a moment) 

-Re{e-^^G'{'tP)xj,X3)L^ = {W^XjWh + i^^\\xXj\\h) +1^ [ XjLxjdx + xT 

-Ay / \tl}\^’^\x 3 \^ dx + aX ( (i/>Xj)(/3Im (i/'x)') - 7^6 

where we have also used the fact that ||x(t)||L2 = 1. Next, we estimate the term 
proportional to H as we did in the proof of Lemma 14.21 and we also use that fact 
that 

f P^e{tpXj)l3^x[i{'ipx-)dx ^\P\ [ \tp\'^\xj\'^ dx. 

In summary, this yields 

-Re(e"®’’G'(V')Xj,Xi)L2 < “ 2 (II^Ailli2 + w2||aiXilli2) + 7M 

-A(7-a|/3|)/ 


m ^2 ^ 

-^Re(e"®’’G'(i/:)xj,Xi)L2 ^ -7(1- 2) '^i^oXj, Xj) l'^ +jp.m 


Thus, 
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in view of definition (12.21) . To further estimate the right hand side of (16.41) , we use 
Holder’s inequality and Gagliardo-Nirenberg to obtain 



We will use the generalized Sobolev-Lieb-Thirring inequality (see [52]) that reads 



provided max{l, d/2} < p < l+d/2. Here ci = ci{d,p) > 0 some absolute constant. 
Choosing p = 1 + i (which requires a > 2/d), we obtain 



Young’s inequality then implies that for any e > 0, there exists a C 2 = C 2 {d, a) > 0, 
such that 


2(t + 1) 


2a 

^ 2 <t +2 


Eii^x 






t=i 




^2 II / ||2<tq: 


L2“+2 + £E(^0XjlXj)L^; 


t = l 


where a = 2 a+ 2 -da ’ ^ “ 2 a+ 2 -da ' Note that both of these exponents 


2 a-\-2 — da 


are 


positive for a < 2 (d- 2 ) ■ Thus, we an appropriate choice of e, we obtain from (16.4|) . 
that 

m ^2 


t=i 


i=i 

C3(A|^|)1+“ 


/yO 




2 aoL 

]/^2(7 + 2 . 


Now, using the estimate from Lemma 16.91 above, we have 


02 

-ER'e(e"*’^G'(V')xi,Xj)L2 < - ^ 


.2»(i-5)-"mI 


i=i 


2(7 Q! 

L 2 .T + 2 . 


This can be estimated further by 


*’’G'('!/(t))xi,Xi)L2 ^ + p(t), 

where ki is as defined above and 


2(7(5; 

l^2cr+2 , 


with C 4 = 04 ( 0 ;, d) > 0 . 
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Now, for '0(f) = i?(f)'0o G A, we have that 

(5 = limsup sup f ||i?(s) 0 o|lL 2“+2 ds ) < oo, 

*->■00 V'oe^ \ ^ Jo ) 

due to Lemma [4.21 and Proposition 15.31 which imply that for '0(f) G A: 

iiv^wiil 2“+2 <ii0(f)iir“<pr“. 

This consequently yields 

Qm ^ + *t 2 , for all m ^ 1 , 

which finishes the proof. □ 

Remark 6.11. In comparison to many classical results on the dimensions of global 
attractors (e.g., EH), the proof above requires the use of the generalized Lieb- 
Thirring type inequality to control the term proportional to A, see [22] for more 
details. 


We expect that a similar analysis can be done to estimate the box dimension of 
the attractor, cf. m for more details. We finally note that a careful analysis of all 
the involved constants in ki, K 2 shows that for a given, fixed a; > 0 , the fraction 


— ) ^ + 00 , 

Kl 


as |n| 


The estimate on the dimension of A thus becomes larger the larger the rotation 
speed. 


Appendix A. Derivation of the kernel of the linear semi-group 


Our starting point for justifying the formula (12. 6 |) for the kernel of the linear 
semigroup, is the following linear Schrbdinger equation 

idtu = Hqu, u\t=o = uo{x). 

where, as before, iJo = + kP)- For this equation, Mehler’s formula yields 

an explicit representation of (the kernel of) the associated semi-group [5]. More 
precisely. 


u { t , x )= j So { t , x , y ) uo { y ) dy , 


where 

(A.l) 

and 


So { t , x , y ) = { 2 Trigo { t )) = exp 

. , sin(a;t) 
go { t ) = 


9 o(t) 


ho(t) 


[ x ^ + y ^)- x-y 


ho { t ) = cos { ujt ). 


It was shown in [4] , that a simple change of variables allows to obtain an analogous 
formula for solutions of the linear Schrbdinger equation with non-vanishing rotation. 
To this end, we write 

xi = cos(Ot)ii -I- sin(nt)i2, X2 = cos { flt ) x2 — sin(nt)ii, 

and xs (if applicable) is left unchanged. Note that this transformation is volume¬ 
preserving and hence does not affect the pre-factor { 2 Trigo{t))~‘^A which ensures 
that 


R<*xR‘^ 


So { t , X , y ) dx dy = 1 , for alH > 0 . 
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Defining the new unknown u(t,x) = u{t,x) a straightforward calculation shows 
that u solves 

(A.2) idtu = Hqu — flLu, u\f^Q = uo{x), 

Substituting the new coordinates into Mehler’s formula (IA.1|) yields 


Snit,x,y) = {2Trigo{t)) 2 exp 


9oit) 


ho{t) 


(x^ + y^) — {cos{flt)xi +sin(f2t)a:2)2/i 


— {cos{flt)x 2 


sm{nt)xi)y 



Here we drop the term x^y^ for notational convenience, since it is unchanged by 
the change of coordinates. In order to finally obtain SQ(t,x,y), i.e., the kernel for 
the dissipative semi-group 5*0(t) associated to (11.71) we replace 1 1->- —in the 
above kernel. In other words, 

Sn{t,x,y) = Sn {-ie~‘^^t,x,y) , 

which, after some algebra, yields (12.61) . In there, the pre-factor in front of the 
exponent is understood in terms of the principal value of the complex logarithm 
via (a -I- ib)'* = and is differentiable for small enough f > 0. 

Next, in order to study the regularizing properties of 5'o(t) for short times, we 
first note that the phase function $ in (12.61) can be decomposed into its real and 
imaginary part, denoted by $ = $1 -I- *<I’2. For Gaussian, i.e. heat kernel type, 
regularity properties of the semi-group 5'o(t) for small f > 0, we require (at least) 
an inverse quadratic decay of the real part $1. To this end, let 

n 

a = tu! cos'd, b = tujsm'd, vj=—. 

OJ 

With this choice of notation we have 


1 




\ sinh(a -I- ib) ) 
Im 


cosh(ci7(a -I- ib)) 
sinh(a -I- ib) 


( sinh(ri7(a -|- ib)) \ , . 

( Sinl.(a + .t) 


[xiyi +X2y2) 


Let us investigate the behavior of the real part $1 of the exponent near t Ri 0. 
Standard trigonometric identities yield 


^ / cosh(^(a + ii,)) 
\ smh(a -I- lb) 


cosh(a7a) cos{wb) sinh(a) cos(6) -I- sinh(n7a) sin('n76) cosh(a) sin(6) 
sin^(&)-I-sinh^(a) 


and 


\ sinh(a -I- ib) 

cosh(ti7a) sin(n76) sinh(a) cos(6) — sinh(cc7a) cos{wb) cosh(a) sin(6) 
sin^(6)-I-sinh^(a) 

for all a,b,xu G M. In case tu = 1, the identity 2 sinh(a) cosh(a) = sinh(2a) yields 


(A.5) 


Re 


sinh(2a) 


cosh(a -I- ib) 
sinh(a + ib) J 2 sin(5)2 -|- 2 sinh(a)^ 
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Now, we can Taylor expand the expressions (IA.3|) . (IA.4I) and (IA.5I) around t = 0. 
Recalling a = twcosi?, b = — twsin'd, and ro = £, the denominator of all three 
terms (TOIi-iTAa]) equals 

sin^(6) + sinh^ (a) = sin^ (tuj sin i5) + sinh ^{tuj cos i5) 
Straight-forward expansion yields 

sin^(taj sini?) -|-sinh^(to; cosi?) = cos(2'd) -|- 0 {t^), 


and thus 


: (l — cos(2i?)) -I- 0{t^). 


sin^ {tuj sin i9) + sinh^ (tuj cos d) 

Here and throughout the constant in the 0(-)-notation only depend on i?, w, and 
il. The numerator of (IA.5I) satishes 

sinh(—2taj cosi?) = —2twcosi5— cos^ i?-f 

On the other hand, the terms in the numerator of (IA.3I) satisfy 

cosh(tO cos d) cos(tO sin d) sinh(tw cos i?) cos{toj sin d) 

= tu; cos'd + cosi?(r2^ cos^ i? — -I- sin^ i9) -|- cos^ i9 + O(t^) 


and 


cosh(tO cos I?) cos(tO sin i?) sinh(tw cos i9) cos(tuj sin i?) 

= cos I? sin^ I? -I- O (t^ ). 

Hence, their sum equals 

tiocos'd + cosi?(3H^ -|- 2w^ cos(2i5) — -I- 0(t^). 

Likewise we can easily check 

cosh(fH cos i?) sin(tH sin i?) sinh(ta; cos i?) cos{tuj sin -d) 

= sin d cos d + 0{t‘^) 

= sinh(tH cos d) cos{tH sin -d) cosh(tu; cos d) sin(tw sin d) 

and hence the numerator of (IA.4I) vanishes up to fourth order. Collecting all the 
expansions so far, we obtain 

$i(t, X, y) = ^ (1 — cos(2'd)) + y^){ — 2tu}cosd — cos^ i?) 

-I- 2(a;iyi -I- X 2 y 2 ){ti^ cosd + cosi?(3H^ -I- 2a;^ cos(2i?) — ta^)) 

+ 0((^2+2/V)) 

= +2{xiyi +X2y2){ - 1 + - 30 ^)) 

+ 0{ix^ + y^)t^)y 

This expression can be further simplified by collecting coefficients of x — y and x + y 
to obtain 

$i(t,x,y) = -^^(^{x-yf{2+^{uj^ + 3n^)t^) + 

{x + y)2i(a;2 - + 0{{x^ + y2)t3)j . 
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Note that cos(i9) > 0 if i9 G (—§, f) and we assume that uj > ^ 0. The term 

0{{x^ +y^)t^) can thus be absorbed in the other coefficients for small enough t. In 
particular, if t < d for some small enough d > 0, the real part satisfies 

^i{t,x,y) ^ , c> 0. 

ct 

Hence the semi-group has the same decay as the heat kernel and indeed satis¬ 
fies (12.71) . Since the exponent F is quadratic in x and y, the derivative of S w.r.t. 
X yields only an extra linear factor. In summary, we find that for small t > 0, the 
absolute value of the kernel is bounded by 

|5'n(t,a:,y)| ^ {t + 0{f))~^ exp - 2/)^(2 + i(a;^ -k 30^)^^) 

-k (x -k - Xsys + 0((x^ -k y‘^)t^)j ^ ■ 
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